Dynamic crossover in the persistence probability of manifolds at criticality 
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We investigate the persistence properties of critical d-dimensional systems relaxing from an initial 
state with non- vanishing order parameter {e.g., the magnetization in the Ising model), focusing on 
the dynamics of the global order parameter of a d'-dimensional manifold. The persistence probability 
Pc{t) shows three distinct long-time decays depending on the value of the parameter £ = (D—2+r])/z 
which also controls the relaxation of the persistence probability in the case of a disordered initial 
state (vanishing order parameter) as a function of the codimension D = d — d' and of the critical 
exponents r\ and z. We find that the asymptotic behavior of P c (t) is exponential for £ > 1, stretched 
exponential for < £ < 1, and algebraic for £ < 0. Whereas the exponential and stretched 
exponential relaxations are not affected by the initial value of the order parameter, we predict and 
observe a crossover between two different power-law decays when the algebraic relaxation occurs, 
as in the case d! = d of the global order parameter. We confirm via Monte Carlo simulations 
our analytical predictions by studying the magnetization of a line and of a plane of the two- and 
three-dimensional Ising model, respectively, with Glauber dynamics. The measured exponents of 
the ultimate algebraic decays are in a rather good agreement with our analytical predictions for the 
Ising universality class. In spite of this agreement, the expected scaling behavior of the persistence 
probability as a function of time and of the initial value of the order parameter remains problematic. 
In this context, the non-equilibrium dynamics of the 0(n) model in the limit n — > oo and its subtle 
connection with the spherical model is also discussed in detail. In particular, we show that the 
correlation functions of the components of the order parameter which are respectively parallel and 
transverse to its average value within the 0(n — > oo) model correspond to the correlation functions 
of the local and global order parameter of the spherical model. 
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I. INTRODUCTION 



After decades of research, understanding the statistic of first-passage times for non-Markovian stochastic processes 
remains a challenging issue. Of particular interest in this context is the persistence probability P c (t), which, for a 
stochastic process X(t > 0) € K of zero mean (X(t)) = 0, is defined as the probability that X does not change 
sign within the time interval [0, t]. In terms of P c (t), the probability density of the first time at which the process 
crosses X = (zero crossing) is — dP c (t)/dt. This kind of first passage problems have been widely studied by 
mathematicians since the early sixties [l|-|3(, often inspired by engineering applications. During the last fifteen years 
these problems have received a considerable attention in the context of non-equilibrium statistical mechanics of 
spatially extended systems, both theoretically 0; @ and experimentally @. In various relevant physical situations, 
ranging from coarsening dynamics to fluctuating interfaces and polymer chains, P c {t) turns out to decay algebraically 
at large times, P c (t) ~ t~ 9c , where 6 C is a non trivial exponent, the prediction of which becomes particularly challenging 
for non-Markovian processes @ • 

In statistical physics, persistence properties were first studied for the coarsening dynamics of ferromagnetic spin 
models evolving at zero temperature T = from a random initial conditions Q . In this case the local magnetization, 
i.e., the value ±1 of a single spin, is the physically relevant stochastic process X(t) and the corresponding persistence 
probability turns out to decay algebraically at long times. By contrast, at any non-vanishing temperature T > 0, the 
spins fluctuate very rapidly in time due to the coupling to the thermal bath and therefore the persistence probability 
of an individual spin decays exponentially in time. However, it was shown in Ref. Q that the global magnetization, 
i.e., the spatial average of the local magnetization over the entire (large) sample, is characterized by a persistence 
probability P c (t) — referred to as global persistence — which decays algebraically in time P c (t) ~ t~ 6g at temperatures 
T below the critical temperature T c of the model. The case corresponding to T — T c , which we shall focus on here, 
is particularly interesting because 8 g turns out to be a new universal exponent associated to the critical behavior of 
these systems. The global persistence for critical dynamics has since been studied in a variety of instances [8l-[l0| (see 
also Ref. [l2[ for a recent study of the global persistence below T c ). 

In order to understand how the long-time behavior of the persistence probability at fixed T = T C crosses over from 
the exponential of the local magnetization to the algebraic of the global one, Majumdar and Bray [l3[ introduced 
and studied the persistence of the total magnetization of a d'-dimensional sub manifold of a d-dimensional system, 
with < df < d (a similar idea was put forward in Ref. [l(| and further studied in Ref. [Hj]). The two limiting cases 
d! = and df = d correspond, respectively, to the local and the global magnetization and therefore a crossover from an 
exponential (df — 0) to an algebraic (df = d) decay is expected in the persistence probability as df is varied from to 
d. Interestingly enough, it turns out that as a function of d' the persistence probability P c (t) of the manifold displays 
three qualitatively different long-time behaviors, depending on the value of a single parameter £ = (D — 2 + rj)/z, 
where D = d — d' is the co-dimension of the manifold, z the dynamical critical exponent and ry the Fisher exponent 
which characterizes the anomalous algebraic decay of the static two-point spatial correlation function of the spins. As 



where the exponent Oo(d, d') is a new universal exponent which depends on both d' and d, whereas a\ and b\ are non- 
universal constants. Interestingly enough, from the mathematical point of view, the somewhat unexpected stretched 
exponential behavior for < C < 1 emerges as a consequence of a theorem due to Newell and Rosenblatt @ which 
connects the long-time decay of the persistence probability of a stationary process to the long-time decay of its two- 
time correlation function. In passing, we mention that a behavior similar to Eq. ([1]) was also predicted and observed 
numerically for the non-equilibrium correlation function and the global persistence of spins close to a free surface of 
a semi-infinite three dimensional Ising model with Glauber dynamics |14j . 

The results mentioned above — including the crossover in Eq. ([T]) — were obtained for the critical coarsening 
of a system which is initially prepared in a completely disordered state characterized by a vanishing initial value 
mo = of the magnetization or, in general terms, of the order parameter. Due to the absence of symmetry-breaking 
fields, this initial condition implies that the average m(t) of the fluctuating magnetization at time t over the possible 
realizations of the process vanishes. For mo 7^ 0, instead, after a non-universal transient, the average magnetization 
m(t) grows in time as m(t) cx mot 9 ^ for t <C r m cx m whereas, for t ^> r m , m(t) decays algebraically to zero as 
m(t) oc t~P/^ vz \ These different time dependences are characterized by the universal exponents 9[ s > (the so-called 



a function of ( one finds [l3[ 




for C < 0, 
for < C < 1, 
for C > 1 , 



(1) 



initial-slip exponent [lj|) and 



K = e[ s + p/(vz), 



(2) 
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where /3, v and z are the usual static and dynamic (equilibrium) critical exponents, respectively. In Ref. fl6| . we have 
demonstrated that a non- vanishing value of tuq results in a temporal crossover in the persistence probability P c (t) of 
the thermal fluctuations Sm(t) of the fluctuating magnetization around its average value m(t), such that 

p ,,s jt- 00 for t micI <t<r m , 
\t°°° lor t > r m , 

where f m i c r i s a microscopic time scale. On the basis of a renormalization-group analysis of the dynamics of the order 
parameter up to the first order in the dimensional expansion and of Monte Carlo simulations of the two-dimensional 
Ising model with Glauber dynamics, we concluded that the two exponents 8 and 9^ arc indeed different, with 

&o < Ooo- 

In the present study, we investigate both analytically and numerically the interplay between the crossovers described 
by Eqs. (JXJ) and ([3]), focusing primarily on the case of the Ising model universality class with relaxational dynamics. We 
argue that for £ > the qualitative behavior of the persistence probability in Eq. (fTJ) is not affected by a non- vanishing 
initial magnetization. For £ < 0, instead, we predict a crossover between two distinct algebraic decays characterized 
by different exponents 9 (d,d r ) and 9 00 (d, d'), as in Eq. ©. The conclusions of our analysis of the time dependence 
of the persistence probability are summarized in Tab. HI 







T m < t 


C<o 


Pc(t) ~ £-"o(cW) 


Pc{t) ~ t- 9 ^ d ' d,) 


o < c < i 


P c (t) ~exp(-ai£ c ) 


C> i 


P c (t) ~ cxp(-M) 



TABLE I. Time dependence of the persistence probability P c (i) of a d'-dimensional submanifold of a d-dimensional system 
relaxing from an initial state with non- vanishing order parameter mo, as a function of the parameter £ = (d — d' — 2 + r])/z. 
(See the main text for details.) 



The rest of the paper is organized as follows: In Sec. UH we describe the continuous model we shall study and we 
present a scaling analysis which yields the behaviors mentioned above for the persistence probability P c (t ^ r m ) 
at criticality, see Tab. Q] In Sec. IIII1 we present an analytic approach to the calculation of 9 00 (d,d / ) for the Ising 
universality class with relaxational dynamics. In particular, in Sec. IIII A[ we focus on the Gaussian approximation, 
whereas in Sec. IIII B[ we discuss the effects of non-Gaussian fluctuations. In Sec. IIV1 we present the results of our 
numerical simulations of the Ising model with Glauber dynamics and we compare them with the analytical predictions 
of Sec. Mil In Sec. [V] we study the persistence of manifolds for the 0(n) model in the limit n — > oo and, in passing, 
we discuss the connection between the non-equilibrium dynamics of this model and of the spherical model. Our 
conclusions and perspectives are then presented in Sec. IVII 



II. MODEL AND SCALING ANALYSIS 



Here we focus primarily on the Ising model on a d-dimensional hypercubic lattice, evolving with Glauber dynamics 
at its critical point and we study the persistence properties of the associated order parameter both analytically and 
numerically. The universal aspects of the relaxation of this model are captured by the case n = 1 of so-called Model 
A [l7j for the n-component fluctuating local order parameter </? Q (x, t), a = 1, . . . ,n (e.g., the coarse-grained density 
of spins at point x = (xi, . . . , xj) in the Ising model): 



d t (f(x,t) = - 



6H[<p] 



?7(x,t), 



(4) 



where 77 (x, t) is a Gaussian white noise with (T) a (x.,t)) = and (r? a (x, ^(x', t')) = 2T5 a p5 d (x - x')6(t - t') and T 
is the temperature of the thermal bath. In what follows we consider the critical case T = T c . In Eq. (TJ}, the friction 
coefficient has been set to 1 and % is the 0(n)-symmetric Landau-Ginzburg functional: 



H[ V ] = J d d x 



^(V^ + ^ + f^ 2 ) 2 



(•5) 
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where J d d x = J J\i=i d-Xi, tq is a parameter which has to be tuned to a critical value ro. c in order to approach the 
critical point at T — T c (here r 0iC = 0), and go > is the bare coupling constant. 

At the initial time t = 0, the system is assumed to be prepared in a random spin configuration with a non-vanishing 
average order parameter Mq along the direction in the internal space which will be denoted by 1, [ip a (x, 0)]o = MoS a i 
(a = 1, . . . ,n) and short-range correlations [<p a (x, 0i)<pp(x!, 0)]o — M^SaiSpi = r " 1 (5 C[( g(5' i (x — x'), where [. . .]o stands 
for the average over the distribution of the initial configuration. It turns out that Tq -1 is irrelevant in determining 
the leading scaling properties [lj| and therefore we set t^ 1 =0. Here we first focus on the case n = 1 of the Ising 
universality class, whereas in Sec. |V]we shall extend our analysis to n > 1. The stochastic process we are interested 
in is the dynamics of the fluctuations of the magnetization around its average value, i.e., 

ip(x,t) = ip{pc,t)-M{t), where M(t) = (<p(x, t)) (6) 

is the average magnetization and (. . .) stands for the average over the possible realizations of the thermal noise rj. 
Assuming the total system to be on a hypercubic lattice of large volume L d , the total (fluctuating) magnetization of 
the d'-dimensional sub manifold is given by 



M(x d >. 



1 f d ' 

■ > x d>t) = -j-jT YidXil(>(xi,...,Xd',Xd> + i,...,Xd,t) , (7) 



and it depends on the remaining D = d—d! spatial coordinates, which shall be denoted by the vector r = (n, . . . ,ru) = 
(Xd' + l, ■ ■ ■ ,x d ). 

The analytical calculation of the persistence probability P c (t) of the process M(r,t) relies on the observation that 
the field M(r,t) is a Gaussian random variable for d! > 0. Indeed, for a d'-dimensional manifold of linear size 
L, M(r,t) is the sum of L d (d! > 0) local fluctuating degrees of freedom, e.g., spins in a ferromagnet, which are 
correlated in space only across a finite, time-dependent and growing correlation length £(t). In the thermodynamic 
limit L > ^) the number ~ [L/^(t)] d of effectively independent variables contributing to the sum becomes large and 
the central limit theorem implies that M(r,t) is a Gaussian process, for which powerful tools have been developed 
in order to determine the persistence exponent 0, [l8j]. In particular, the Gaussian nature of the process implies that 
P c (t) is solely determined by the two-time correlation function Cj^{t,t') = (M(r,t)M(r,t')). For t' < t <C r m , the 
correlation function C^{t, t') coincides with the one of the case with vanishing initial magnetization r m — > oo, which 
was studied in Ref. [13j |. The corresponding persistence probability P c (t) displays the behavior summarized in Tab.|T] 
for t <C T m . However, a non- vanishing average Mq of the initial order parameter affects the behavior of C^(t,t') as 
soon as t',t ~ r TO . In the long-time regime t > t' ^> r m , one can take advantage of the results presented in Ref. [l9| 
for the scaling behavior of the two-time correlation function of the Fourier transform ^(Qi t) = / d d xe jC ^ x t) of 
the local fluctuation ip(x,t) of the magnetization [see Eq. ©]. Here and in what follows, the Fourier transform of 
^>(x, t) will be simply denoted by using Q = (Qi, . . . , Qd) as an argument of tp and we will also assume t > t' . This 
yields the following scaling form for the Fourier transform M(q, t) oc "0(Q = (0, ■ ■ ■ ,0, qi, ■ ■ ■ , 9d), t) of M(r, t): 

(M(q,t)M(-q,i')) = ^(^) e ' fc (q Z (t - O, |) ' ( 8 ) 

where q = (q u ■ ■ ■ ,q D ), q = |q|, and 9 = -f38/{vz) = -{d+2 - rj)/(2z) [H. With the prefactor (t/t') 6 ' 1 explicitly 
indicated, the function fc is such that fc{v,u — > oo) = fc,oo{v) for fixed v. The correlation function Cj^(t,t') which 
determines the persistence probability is then obtained from the correlation function in Eq. ([S]) by integration over 
momenta 

C M (t,t') = {M(r,t)M(r,t'))= J f{ ^|<M(q,t)M(-q,t')> • (9) 

i—l 

The scaling form (JSj> has the same structure as the one discussed in Ref. [l3| for mo = 0, the only differences being in the 
specific value the exponent 9 and in the form of the scaling function fc{v,u). Accordingly, the argument presented 
in Ref. [HI — which is actually independent of these differences — leads to the conclusion that the persistence 
probability for f>r m behaves as: 



Pc(t) 



t -6 OB {d',d) for £ < o 

exp (-aii c ) for < C < 1, (10) 
exp (—bit) for £ > 1, 
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where the non-universal constants a\ and b\ are the same as in Eq. (p}. Indeed, in the case £ oc D — 2 + 7/ > 0, 
the integral in Eq. © for the equal-time correlation function Cj^(t — t',t') oc J d D q q~ 2+v would diverge at large 
wavevectors in the absence of a large-wavevector cutoff A. Accordingly, for times large enough compared to the 
temporal scale t\ ~ A _z associated to that cutoff the integral is dominated by the quasi-equilibrium regime t ~ t' 
(i.e., t, t' 3> t\ with arbitrary r = t — t'), which is actually independent of the initial value mo of the magnetization 
and is therefore identical to the case £ > discussed in Ref. [l3| and corresponding to mo = 0. (Note that, strictly 
speaking, the application of the theorem by Newell and Rosenblatt 0] to the case < £ < 1 actually provides the 
corresponding expressions in Eqs. (fTU|) and ([T]) as upper bounds, see Ref. H3 for details.) 

The exponent 0oo(d', d) which describes the algebraic decay of P c (t) for C < 0, instead, is a new universal exponent 
and it is expected to differ from 9o(d' , d) in Eq. |T]). Indeed, the resulting persistence exponent depends on the specific 
value of 9 and on the specific form of the scaling function in Eq. l[8]l. which both depend on having t, t' <C T m (mo = 0) 
or t, t' > T m (mo ^ 0). 

Accordingly, we shall focus below on this novel algebraic behavior emerging for £ < 0, which is characterized by 
the exponent 9 oc (d,d') (see Tab.QJ. 



III. ANALYTIC CALCULATION OF THE PERSISTENCE EXPONENT 

In the case we are presently interested in, D — 2 + rj oc £ < and the integral in Eq. © — necessarily defined 
with a large wave- vector cutoff A — is convergent for A — > oo so that C^(t,t) = (M 2 (r,t)} [and therefore Cj^(t,t')] 
has a well-defined A-independent limit which can be taken from the outset, provided that the expression below are 
understood to be valid for times much larger than the time scale t\ [— t nlicr of Eqs. (fTJ) and ([3])] associated to that 
cutoff. 

In order to study the persistence properties it is convenient to focus on the normalized process X(r,t) = 
M(Y 1 t)/(M 2 (Y,t)) 1 / 2 associated to M(r,t) (see, for instance, Ref. 0), which is characterized by a unit variance. The 
scaling form ([5]), together with Eq. ©, implies that the correlation function of X(r,t) has the following scaling form 
for t > t' > r m 



(X(r,t)X(r,t'))=^-j Foo(*/f), (H) 

where Foq^x) is a non-constant function, with the asymptotic behavior F QO (x — > oo) ~ const ^ and 

fx 00 (d,d') = -(e-i) + ^ (i2) 

d + D 

= 1 + —- ^ 

Even though the Gaussian process described by Eq. (|II[) is not stationary, it becomes such when the logarith- 
mic time T = log* is introduced. In fact, this yields (X(T)X{T')) = e-'^ T - T V 00 (e T - T ') = C st (T - V) with 
C s t (AT) ~ exp[— fJ,oo(d, d')AT] for large AT, which implies Q that the persistence probability P c decays exponen- 
tially in logarithmic time, P C (T) ~ exp[—9 00 (d,d')T}, i.e., an algebraic decay P c {t) ~ t-6<*>(d,d) m ^ erms f 
original time variable. In addition, if F^x) is actually independent of x, i.e., if it is a constant, the process can be 
mapped onto a Markovian process for which 9 oa (d,d') = /J<oo(d, d'). In a sense, /ioo provides a sort of "Markovian 
approximation" of the exponent 9^. In the generic case, however, the exponent 9 OQ (d,d') depends in a non-trivial 
fashion on both the exponent [ico(d, d') [see Eq. (|lip] and the full scaling function F OQ (x), which, for the present case 
D 0, is known only within the Gaussian approximation of Eqs. (fj| and ([5]) [19] (whereas for D = 0, an expres- 
sion for the first correction in the e-expansion about the spatial dimension d = 4 — e is available [l9l]). It actually 
turns out (see further below) that the function F oc (x) is non-trivial already within the Gaussian approximation, a 
fact that makes the calculation of the exponent 9oo(d,d') a rather difficult task. Here we present a perturbative 
expansion of this exponent which is valid for small co-dimension D « 1. Indeed, one notes that F x (x) reduces 
to a constant for D = and therefore the persistence exponent 9 oc (d,d') is given by its Markovian approximation 
oo (d,d') = fj,oo [1, [l6|. For small D one can take advantage of the perturbative formula derived in Refs. 0, Q in 
order to expand the non-Markovian process for D ^ and its persistence exponent 9^ around the Markovian one 
corresponding to D = 0. 
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A. The Gaussian approximation 



The Gaussian approximation discussed below is actually exact in dimension d > 4 and it is obtained by neglecting 
non-linear contributions in ip to the Langevin equation Q once it has been expressed in terms of ip[x,t") and m 2 = 
goM 2 /2 (see, e.g., Ref. [19(). This yields the evolution equation 

<9 t V(x, t) = [V 2 - m 2 (£)]V(x, t) + ??(x, t) 

1 (14) 
where d t m(t) + —m (t) — , 
o 

with the initial condition m(t = 0) = mo. In order to determine the persistence exponent, we need to calculate the 
correlation function C-^ of the order parameter of the manifold [see Eqs. (J7J and (|SJ)]. In turn, Cj^ can be inferred from 
the two-time correlation function of the Fourier components V*(Qi = / d d xe^ x- 0(x, t), where Q = (Qi, . . . , Qd), 
which was calculated in Ref. (see Eq. (59) therein): 

2 e -Q 2 (*+*') /■*< 2 

C Q (t,f) = (V(Q,*)V'(-Q,*)> = F77- w. /-L , 13/2 / dti(t! +r m ) 3 e 2 « 41 , (15) 

where i< = min {t,f} and r m = 3/(2m§). The correlation function C~(£,7j') = (M(r,i)M(r,i')) in space follows from 
the integration (0) of the correlation of M(q, t) cx V>(Q = (0, ... ,0, (ft, ... , qo),t) [see Eq. (J7J]: 

C^(M') =<M(r,i)M(r,0) = / J] ^|(7 Q=(0jq) (*,*') 

(16) 



2cdt, 



l-D/2 ~t 



[(£+ + l)] 3 / 2 7o 



d^ (t 1 + lf(t + t' -2h)- D/2 



where, on the first line, we used the notation Q = (0, q) = (0, • • • , 0, q±, Q2, • • • , Qd) and, on the second, we introduced 
cd = (47r) _£) / 2 and the dimensionless time variables t = t/r m and t' = f'/r m , assuming i' < t. The correlation 
function (X(r,t > t')X(r,t')) of the normalized process X(r,t) = M{r,t)/(M 2 (r,t)} 1 / 2 is therefore given by 

(X(r t)X(r t')) = C m(M0 = / 1 \ - ( " +2)/4 fl71 

' ' yJCj^(t, t)Cj^(t' , f) \t'J [/d/ 2 (M)/d/ 2 (M')] 1/2 ' 

where 

I a (x,u)= [ dv(l + uv) 3 [l +x(l - 2v)]- a . (18) 



Note that in order for C^(t,t) to be defined, Id/2(1,w) ~ J 1 dv(l — v)~ D / 2 has to be convergent, which requires 
D < 2, consistently with the assumption £ = (Z? — 2)/2 < (77 = and z = 2 within the present approximation). In 
contrast to the case D = of the global persistence, discussed in Ref. [l6| [see Eq. (6) therein] , there is no suitable 
choice of a function L(t) such that the correlation function (fTTf for D ^ takes the form of a ratio L(t')/ L(t), which 
corresponds to a Markovian process. Accordingly, the Gaussian process {X(t)}t>o f° r D 7^ displays significant 
non-Markovian features even within the Gaussian approximation. 

In the asymptotic regime t > t' >■ 1 we are interested in I a (x, u 3> 1) ~ u 3 dvv 3 [l + x(l — 2v)]~ a and therefore 
the correlation function (fT"7|) takes the form 

(X{r, t)X(r, t')) ~ (t/<')" (2+D/4) -4(t'/t) , for <> t' » r m , (19) 

where 



A(x)=A D f dvv 3 [l + x(l-2v)]- D/2 , (20) 
Jo 

and A D = [2~ D / 2 x 6r(l - L>/2)/T(5 - D/2)}- 1 = (8 - D)(6 - D)(4 - D)(2 - £>)/(3 x 2 5 --°/ 2 ), with finite .4(0) = 
and, by definition, 4(1) = 1. The correlation function (|19l) has the form (jlll) with d — A, z — 2, i.e., 
/ioo = 2 + D/4 and Foo(t/t') — 4(t'/t). The scaling function A(t' /t) is indeed a non-trivial function of its argument 
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and it reduces to a constant only for D = 0, i.e., for vanishing codimension. As anticipated, one can take advantage 
of this fact in order to determine oo (d,d / ) perturbatively, by expanding around the Markovian Gaussian process for 
D = according to Ref. 9]. First, one introduces the logarithmic time T = logi and expands (X(r,T)X(r,T')) for 
small D: 



with fi x = 2, 



and 



Ax(x) 



(X(r,T)X(r,T')) = c -^(t-t') ^-(t-t')) ? 



A(x) = x D/4 A(x) = 1 + DAx(x) + 0{D 2 ) , 



~ + | + y +V 3 ~ ^log(2y- 1) + (y 4 - l)log(j/ - 1) - y 4 logy 



(21) 
(22) 

, where y= (x^ 1 + l)/2. (23) 



The function A\{x) is responsible, at the lowest order in D, for the non-Markovian corrections to O 
can be calculated using the perturbation theory of Ref. @ : 



TZ 



= 1- D 



2^c 



Ax- — 

o (1 



>)3/2 



0{D 2 ). 



IZHoo which 



(24) 



With the explicit expression for A\(x) given in Eq. (|23[) . a straightforward numerical integration yields 

TZ = l + D x 0.848- ■■ + 0{D 2 ), (25) 



which renders 



6W(d,d' = d-D) = 2[l + D x 0. 



0(D 2 )} 



(26) 



as a Gaussian estimate of 0oc. In particular, for D = 1 
the value 



the case numerically studied in Sec. IIVI below — it takes 



0W(d,d' = d-l) = 3.69.. 



(27) 



Alternatively, a different (yet equivalent) numerical estimate of 0^ can be obtained by expanding up to first order in 
D only A(x) = 1 + -D^li(x) + 0(D 2 ) on the rhs of Eq. (fl9| while keeping the full D-dependence of the value of fioo = 
2 + D/4, which corresponds to the Markovian approximation. Accordingly, one has 9 00 (d,d') = TZ = (2 + D/A)1Z 
where the ratio TZ is given here by Eq. in which Tioo i-> Moo = 2 + O(-D) and ^4i(a;) H> -4i(:r) = ^i(x) — (lnx)/4. 
With these substitutions one finds TZ = TZ — D/8 + 0(D 2 ) and therefore 



9M(d,d' =d-D) = 2(1 +D/8)[1 + D x 0.723 • ■■ + 0{D 2 )}, 



(28) 



which, as expected, has the same series expansion as 9^ (d, d' = d — D) up to the first order in D and provides the 
numerical estimate 



9^(d,d' =d- 1) =3.87. 



(29) 



for D = 1. These expressions for the persistence exponent 9 oc (d,d — D) have been derived accounting only for 
the effects of Gaussian fluctuations of the order parameter and therefore they are increasingly accurate as the non- 
Gaussian fluctuations become less relevant, i.e., as the spatial dimensionality d of the model approaches and exceeds 
4. Accordingly, the numerical estimates #lo 2 ' (d, d — D) in Eqs. (|2"B"|) and (|28p are expected to be increasingly accurate 
as d increases for a fixed small codimension D. In the next section we shall compare these analytical predictions, 
extrapolated to D = 1 and provided by Eqs. (f27|) and ([29| . to the results of Monte Carlo simulations of the Ising 
model with Glauber dynamics in spatial dimensionality d = 2 and 3, discussed in Sec. IIVI Figure Qla) summarizes 
the available estimates of 9aa(d,d — D) as a function of the codimension D and of the space dimensionality d of the 

model. The solid and the dashed lines correspond to the estimates 0^ (d, d—D) [Eq. (f2"o) ] and 0oo (d, d—D) [Eq. (f2"5)) ] , 
respectively, derived within the Gaussian model. The vertical bars for D = and D = 1 indicate the corresponding 
Monte Carlo estimates (see Ref. [r| and Sec. IIVI below, respectively) for the Ising model with Glauber dynamics 
in d = 2 (grey) and d — 3 (black). As expected, the Gaussian approximation provides a rather accurate estimate 
of the actual value of 0^, in d = 3. For comparison we report here also the codimension expansion of 0o(d,d') - 
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FIG. 1. Persistence exponents 6oa,o(d,d — D) of the global order parameter of a manifold of codimension D as a function of 
D, in the case of (a) non-vanishing and (b) vanishing initial value of the order parameter. The solid and the dashed curves 
in the two panels correspond to the estimates 9^ [Eqs. (|26p and (|33f) ] and 6^ [Eqs. (|28[) an d l|34[) ]. respectively, derived 
within the Gaussian model universality class with relaxational dynamics [Eq. Q]. For D > 2, i.e., > (shaded areas) the 
asymptotic decay of the persistence probability of this model is no longer algebraic. The vertical bars for D = (see Ref. fl6j ] 
and Refs. 20, 21] for (a) and (b), respectively) and D = 1 (present work) indicate the corresponding Monte Carlo estimates of 
the persistence exponents for the Ising model universality class with Glauber dynamics in d — 3 (black) and d = 2 (grey). On 
panel (b), the black and grey circles for D = 1 indicate the preliminary Monte Carlo estimates of Ref. 



corresponding to the case mo = — which can be determined as detailed above for 9^ on the basis of Eqs. (fT7| and 
(fT8f in the limit i~' < i < 1. In this case, I a (x,u< 1) ^ ft du[l + x(l - 2t>)]~ a = [(1 + x) 1 " - (1 - x) 1 - a ]/[2(l - a)x] 
and therefore 

(X(r,t)X(r,t')) - (t/t')- {2+D)/4 B{t'/t), for r m » t > t' , (30) 

where 

(l + sQi-^-q-a:) 1 -^ 

= 2 !-^/ 2 X ' ( } 

with finite B(0) = 2 D / 2 (1 - D/2) and, by its very definition, B(l) = 1. The Markovian approximation to 9q is given 
by hq = 1/2 + D/4 and corresponds to B(a;) = 1 in Eq. fj30[) . i.e., to having D = 0. This expression for fio agrees with 
the hyperscaling relation Q 

Mo = -(^-l) + C/2, (32) 



which can be derived by the same arguments which lead to Eq. (fT2|) , taking into account that in this case mo = the 
correlation function of the order parameter fluctuations scales as in Eq. ((5J) with 9 replaced by the initial-slip exponent 
6>i S = 6[ s — (2 — 7] — z)/z For the present Gaussian model 6*i s = and r/ — 0. (Note that with the identification 
X c /z i— > — (9i s — 1) + C, Eq. (|3^|) renders the expression for the persistence exponent no which is reported after Eq. (8) 
in Ref. [13[ in terms of A c .) Taking advantage of Eq. (|24[) and the series expansion of B for small D one finds 

9[ 1] {d,d-D) = ^[l + Dx(V2- 1/2) + 0(D 2 )} , (33) 

which coincides with the case d = 4 of the expression derived in Ref. [13| for the 0(n oo) model. An equivalent 
estimate of 9o is obtained as discussed above for #oo : 

9[ 2] {d, d - D) = i(l + D/2)[l + £> x (V2 - 1) + 0(i? 2 )] . (34) 

Figure [ljb) summarizes the available estimates of 9o(d,d — D) as a function of the codimension D and of the space 
dimensionality d of the model. As in panel (a) of the same figure, the solid and the dashed lines correspond to the 
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estimates 9^ (d, d — D) [Eq. (|33|) ] and # 2 ' (d, d — D) [Eq. (l34l) ]. respectively, derived within the Gaussian model. The 
vertical bars for D = and D = 1 indicate the corresponding Monte Carlo estimates (see Refs. [2(| [2l| and Sec. IIVI 
below, respectively) for the Ising model with Glauber dynamics in d = 2 (grey) and d = 3 (black). On panel (b), the 
gray and black circles for D = 1 indicate the corresponding preliminary Monte Carlo estimates reported in Ref. [l3| . 
which turn out to be, respectively, marginally compatible and significantly different from the ones in d — 2 (grey) and 
d = 3 (black) presented below in Sec. IIVI 



B. Beyond the Gaussian approximation (perturbatively) 

According to the analytical predictions within the Gaussian model indicated by the solid and dashed lines in 
Fig. [I] the persistence exponents 9^ q (d, d — D) increase upon increasing the codimension D > 0, until D reaches 
the value D — 2, above which (shaded areas in Fig. [T]) the relaxation of the persistence probability is no longer 
algebraic. The contributions of non-Gaussian fluctuations to the global persistence exponents and 9q within the 
0(n) universality class have been calculated analytically and perturbatively in Re fs. [Ill and @, 0j, respectively, only 
for the case D = 0. (For D > the persistence exponent 9q is studied in Ref. [13[ only in the limit n — > oo, see 
Sec.lVlbelow.) For the Ising universality class (n = I) these corrections decrease #oo,o compared to its Gaussian value, 
as the spatial dimensionality d of the model decreases below the upper critical dimensionality d = 4. On the basis 
of these behaviors of 9 00 fl(d,d — D) for d > 4 (Gaussian model) as a function of (small) D > and for D = as 
a function of (small) 4 — d, it is difficult to predict analytically the qualitative dependence on D of 9 OQ ^(d,d — D) 
at fixed d resulting from the combined effect of these two competing trends as a function of D and d. The same 
problem arises in the more general case of the 0(n) universality class with generic n. In this respect, it would be 
desirable to extend to the case D ^ the analysis of the contribution of non-Gaussian fluctuations to the persistence 
exponent of the 0(n) universality class, beyond the n — > oo limit, following Ref. 16| and Refs. 0, H| for the cases 



toq 7^ and m = 0, respectively. This requires the knowledge of the analytic expression (e.g., in a dimensional 
expansion around the upper critical dimensionality d = 4) of the wave-vector dependence of the correlation function 
of the order parameter Cq(£, t') beyond the Gaussian approximation [Eq. (|15l) ]. At present, however, such an analytic 
expression within the 0(n) universality class with Too ^ is available only for Q = (see Refs. [IH, (22[) and, in the 
limit n — > oo, for transverse fluctuations (see Ref. [22]; in Sec. IV Bl we shall comment on the relation between the 
0(n — > oo) model and the spherical model investigated, in this context, in Refs. [23], H3|). For a vanishing initial 
value of the order parameter too = 0, instead, the expression of Cq(£, t') beyond the Gaussian approximation and for 
finite n is known for generic Q only up to the first order in the dimensional expansion around d = 4 (25| and up to 
the second order for Q — (2*1 (see Ref. [27| for a review). However, we point out here that in order to observe a 
non-trivial interplay between the effects of a non- vanishing codimension D ^ and those of non-Gaussian fluctuations 
for d < 4 one would need to account for higher-order terms in the expansion of 1Z around a Markovian process, unless 
the dependence of the correlation function Cci(t,t') on the dimensionality d is known non-perturbatively as in the 
case n — > oo discussed in Ref. [l3[ and in Sec. |V] below. Indeed, the general structure of the correlation function 
of the normalized process is (X(t)X(t')) = (t/t')-^ d ' D ^F(t' /t;d,D) with t > t', where the exponent n(d, D) and 
the function F depend on the specific process under study and fi(d,D) is determined such that F(0;d,D) is finite 
and non-zero [note that, by definition, F(l;d,D) = 1]. If F(x;d,D) turns out to be independent of x [and therefore 
F(x; d, D) = I], then the process is Markovian with persistence exponent /x(<i, D) which, in turn, is related to known 
exponents via hyperscaling relations such as Eqs. p^|) and (|3"2")l . In the limit n — > oo this is actually the case for 
D = and generic d (c.f. Sec. EJ, i.e., F(x;d,D = 0) = 1. Accordingly, the expansion for small D of F[x\d,D) 
takes the form F(x;d,D) = 1 + DdDF(x;d,0) + (D(D 2 ) and one can use the formula presented in Refs. 0, Q [see 
Eq. (IM1) above] in order to calculate the correction 1Z(d,D) = I + Did + 0(D 2 ) to [i(d,D) which determines the 
persistence exponent 9(d,D) = /j,(d,D) x lZ(d,D). Here the deviation from the Markovian evolution is controlled in 
terms of the (small) parameter D and 1Z has a non-trivial dependence on the dimensionality d of the model. This 
approach was adopted in Ref. [l3[ for to = and, c.f., in Sec. [V]of the present work for m 7^ 0. In the case of 
finite n, instead, the correlation function Cq (£,£') for d < 4 is typically known in a dimensional expansion around 
d = 4, up to a certain order in e = 4 — d. As a consequence, F(x;d,D) = F(x;A,D) — eddF(x;A, D) + 0(e 2 ), 
where the lowest-order term F(x;4,D) corresponds to the Gaussian approximation, discussed above for n = 1. 
[The line of argument presented below actually applies also to those cases in which the first correction term to 
Cci(t,t') and therefore to F(x;d,D — 0) is of 0(e 2 ) [7|, [9J.] In the generic case, F(x;4, D) as a function of x is not 
constant, resulting in a non-Markovian process even for e = 0. However, the process turns out to be Markovian 
for D = 0, i.e., F(x;d = 4, D = 0) = 1, and therefore one is naturally led to perform a codimension expansion of 
F(x;4:,D) = 1 + Dd D F(x;4,0) + 0(D 2 ) and, for consistency, of d d F(x;±,D) = d d F{x;A,Q) + 0{D). Accordingly, 
one has F(x; d, D) = 1 + D8dF(x; 4, 0) — eddF(x; 4, 0) +0(D 2 , e 2 , eD), where the deviations from the non-Markovian 
evolution are jointly controlled by D and e. Due to the fact that the perturbative expression for 1Z [see, e.g., Eq. (|24l) ] 
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d = 2 


d = 3 


x> ( n s) 
n, 

77 (ng) 


1.03+0.414D 
1.05 + 0.7237) 


1.01 + 0.4147) 
1.03 + 0.72373 


»7 

z 


1/4 
2.17 
0.38 


0.03 
2.04 
0.14 


Mo 

[-Loo 


0.216 + 0.237) 
1.46 + 0.237) 


0.377 + 0.2457) 
1.74 +0.2457) 



TABLE II. Estimates of the universal ratios ftooo [up to 0(7) 2 ), see Eq. (3SJ], of the critical exponents (see, e.g., Refs. [271, |28|] ) . 
and of the resulting Markovian approximations /xo.oo [Eqs. (|32[) and ()13|) ] for the persistence exponents #o,oo within the Ising 
universality class in two and three spatial dimensions. The analytical estimates 0qoo = ^Ooc x Mo.oo of #o,oo agree with those 
of Refs. [1, [3 for D = 0, whereas for 7) = 1 they yield the values reported in, c.f., Tab. Hill 



is valid up to the first order in the deviation from the Markovian evolution, all the terms of 0(D 2 , e 2 , eD) in F(x; d, D) 
can be neglected when calculating 

K{d = 4 - e, D) = 1 + D A + e B + 0{e 2 ,eD, D 2 ), (35) 

where the coefficients A and B are given by the integrals of 9o-F(a:;4, 0) and — dd,F(x; 4, 0), respectively, according 
to the rhs of Eq. ([2"4")l in which, for consistency, the value fj,(d = 4, D = 0) of fj, at the lowest order in e and D 
enters. Due to the double expansion in Eq. (|35|) . the lowest-order correction to n(d,D) which results from both a 
finite codimension D and non-Gaussian fluctuations for d < 4 is given by the superposition of the two corresponding 
corrections taken separately, which — for the 0(n) universality class — can be inferred from Ref. [l3[ and the present 
study (D > 0, e = 0), and from Refs. 0,1,13 (D = 0, e > 0), respectively. 

Here we focus on the 0(n = 1), Ising universality class and we consider both the cases tuq = and niQ ^ 0, for 
which the constants A and B in Eq. ((35)) actually take different values. Comparing the expression (|35p of 1Z for e = 
with Eqs. I|28p and (f3"4")l one readily finds that 

|V2- 1 = 0.414... for m = 0, 

I 0.723 ... for m ^ . [ ' 

Analogously, the coefficient B for the cases mo = and mo ^ can be determined by comparison with the results 
for K{A -e,D = 0) presented in Eq. (12) of Ref. Q1] and in Eq. (19) of Ref. d, respectively: 

B= f 0.00754... for m = , 

[0.0273... for m ^ . J 

(Note, however, that for mo = the first term of the expansion of ft (4 - e, D = 0) - 1 is of order e 2 [|.) These values 
result in 



^ (4 _ £ J 'M> =l + 7?xO.414... + e 2 xO.OO754... + 0(73 2 ,e 2 Ae 3 ) for m = , ^ 

- i + x 0.723... + e x 0.0273... +0{D 2 1 eD,e 2 ) for m 7^ , 



which clearly show that the corrections due to a finite codimension are quantitatively more relevant than those due 
to non-Gaussian fluctuations. In particular the dependence of the latter on the dimensionality d — 4 — e is weak 
enough that a simple extrapolation of Eq. (f3"5)) to e = 2 and e = 1 should provide quantitatively reliable estimates 
of ft as a function of D for d — 2 and d — 3, respectively. These estimates, denoted by the superscript (ng), are 
reported in Tab. [Til The Markovian approximation ^{d, D) of the persistence exponents in the two cases mo 7^ and 
m = 0, is respectively given by Eqs. (flU)) and (f3"2"]> . where the critical exponents 77, z, and 9- ls in spatial dimension 
d = 2 and 3 (see, e.g., Refs. [23, [28|]) take the values which are reported in Tab. [TTJ together with the resulting 

expressions for /io,oo as a function of D. The persistence exponents 9q and 9^ can be estimated as ^q" 8 ' = ft ng ' x /i 
and #to 8; ' ) = ftoo 8 '* x /Zoo, which are reported as dash-dotted lines in the two panels of Fig. H for d = 2 (gray) and 
d = 3 (black) together with the Monte Carlo and Gaussian estimates anticipated in the corresponding Fig. Q] The 
comparison among all these estimates will be presented further below in Sec. IIV CI 
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FIG. 2. Persistence exponents 6oa,o(d,d — D) of the global order parameter of a manifold of codimension D as a function of 
D, in the case of (a) non- vanishing and (b) vanishing initial value of the order parameter. The dash-dotted lines in (a) and (b) 
correspond to the analytical estimates of #oo,o which account perturbatively for the effects of non-Gaussian fluctuations (c.f., 
Sec HUB)) within the Ising universality class with Glauber dynamics in d = 3 (black) and d — 2 (grey). These lines terminate 
on the right at D — 2 — r\, as the asymptotic decay of the persistence probability of this model is no longer algebraic for a 
codimension larger than these values. Apart from the addition of these dash-dotted lines, this figure is the same as Fig. [1] 



In order to improve on the analytical estimates of the persistence exponents and to go beyond the linear dependence 
on d and D expressed by Eq. ([38]) (see also Tab. HI|) . one would need first of all an expression of 1Z which accounts for 
non-Markovian corrections beyond the leading order, i.e., an extension to higher orders of the perturbation theory 
developed in Refs. [7|,l9|. Secondly, an analytic expression for the contribution of non-Gaussian fluctuations to Cq(£, if) 
would be required. In the present work, instead of pursuing this strategy, we shall investigate the general features 
of the Independence of the persistence exponent 9^ on the basis of the Monte Carlo results presented in the next 
section and of the analytical study of transverse fluctuations in the 0(n — > oo), presented in Sec. IVl 



IV. MONTE CARLO RESULTS 

In order to test our prediction of a temporal crossover in the critical persistence of manifolds for ( < (see Tab. [TJ 
as well as the theoretical estimate of 9oo as a function of d for D = 1 — see Eq. (125p and Tab. [TT] — we studied via 
Monte Carlo simulations the ferromagnetic Ising model with spins Sj = ±1 and Glauber dynamics in two and three 
spatial dimensions. 



A. Line magnetization within the two-dimensional Ising-Glauber model 

We first focus on the two-dimensional case of a L x L square lattice with periodic boundary conditions, for which 
the size L ranges from 64 to 256 and it is such that the data presented below are not appreciably affected by the 
expected finite-size corrections. We study the temporal evolution of the fluctuating magnetization of a complete line 
(row or column, i.e., d! = 1) selected within the lattice and we calculate the associate persistence probability. Taking 
into account the values of the exponents z and r\ reported in Tab. [ill this choice corresponds to ( = (d— d' — 2 + rj)/z ~ 
—0.345 < and therefore the persistence probability of the line magnetization is expected to decay algebraically at 
large times. The system is initially prepared in a random configuration with N + up and 7V_ down spins, where 
N± = L 2 (l ± mo)/2. Then, at each subsequent time step, a site is randomly chosen and the move Si i— > — is 
accepted or rejected according to Metropolis rates corresponding to the critical temperature T — T c . One time 
unit corresponds to L 2 attempted updates of spins. The determination of the persistence probability P c (t) of the 
fluctuations requires also the knowledge of the global magnetization M[t) = L~ 2 (J2iSi), which we obtained by 
averaging over 2000 realizations of the dynamics. For each of these realizations, we also choose a new random initial 
condition with fixed magnetization mo. The persistence probability is then computed as the probability that the 
fluctuating magnetization of the line L _1 J^ieiinc s i ~ M(t) has not changed sign between t — and time t. This 
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(a) (b) 

FIG. 3. Persistence probability P c {t) of the magnetization of (a) a line and (b) a plane within, respectively, the two- and 
three-dimensional critical Ising model with Glauber dynamics for the two extreme values mo = (top data set) and mo = 1 
(bottom data set) of the initial magnetization mo- These data sets correspond, respectively, to the two different regimes t <r m 
and t ^> r m , within which P c (t) decays algebraically at large times but with distinct exponents O^ 10 ^ (top) and #S 1C ' ) (bottom), 
the best fit values of which are indicated by the associated straight lines. These numerical data provide clear evidence for the 
occurrence of the dynamical crossover discussed in Sec. [TT1 



probability is determined on the basis of more than 10 6 samples. In Fig. Efa), we show the results of our simulations 
corresponding to L — 256 and to two different values of the magnetization tuq. The top curve refers to tuq = 0, 
i.e., T m — > oo, such that the system is always in the regime t <C r m investigated in Ref. T3J. These data are fully 
compatible with a power-law decay with an exponent 9 Q MC \d = 2, d! = 1) = 0.78(5) (see also Tab. IIII|) . which turns 
out to be in a rather good agreement with the value ~ 0.72 reported in Ref. [13j ]. The bottom curve in Fig. [3][a), 
instead, corresponds to mo — 1. In this case, r m <C 1 and after a short initial transient, the system enters the regime 
t 3> T m , within which the persistence probability is expected, according to Sec. [Ill to decay algebraically at large times 
with an exponent 6^ . The numerical data in the figure are compatible with such a power-law decay, characterized by 
an exponent { ™ C \d = 2,d' = 1) ~ 4.2(1) which is significantly larger than the one measured in the case mo = 0. 



B. Plane magnetization within the three-dimensional Ising-Glauber model 

In order to investigate the dependence of the persistence exponent on the space dimensionality d at fixed codimension 
D, we extended the investigation presented above (d — 2, D = 1) to the three-dimensional case, focusing on the 
magnetization of a plane (d = 3, D = 1) within the Ising-Glauber model on a L x L x L cubic lattice with periodic 
boundary conditions. The size L ranges from 16 to 128 and it is such that the data presented below are not appreciably 
affected by the expected finite-size corrections. Taking into account the values of the exponents z and 77 reported in 
Tab.HIl this choice corresponds to £ = (D — 2 + rf)/z ~ —0.48 < and therefore the persistence probability is expected 
to decay algebraically at large times. In Fig. ISJb), we report the results of our simulations on a lattice with L = 128, 
for the two extreme values mo = (top data set) and mo = 1 (bottom data set) of the initial magnetization mo, as we 
did in panel (a) of the same figure for the case of the magnetization of a line within the two-dimensional model. At 
large times one clearly observes that the power of the algebraic decay changes when passing from m = to mo = 1. 
In the former case (top curve) t <C T m and the data are compatible with P c (t) ~ f-s (d=3,d =2) an exponent 

9 MC) (d = 3,d' = 2) ~ 1.10(5). This value is significantly larger than the estimate 6 { MC \d = 3,d' = 2) ~ 0.88 
preliminarily reported in Ref. [13| . However, such an estimate was obtained for rather small system sizes L = 15 
and 31 and therefore it might be biased by finite-size effects. The bottom curve in Fig.[3Jb) corresponds to the case 
t ^ T m and the numerical data for the persistence probability still decay algebraically at large times, but with an 
exponent 6^ C \d = 3, d! = 2) ~ 3.7(1) which is significantly larger than in the case mo = 0. 

Intermediate values < mo < 1 of mo, correspond to finite and non- vanishing r m and therefore one expects 
P c (t) to display the two different power-law behaviors described separately above within the two consecutive time 
ranges t <C r m and t ^> r m . Indeed, this is clearly displayed in Fig. |4ja) where we report the time dependence of 
the persistence probability P c (t) of the magnetization of a plane in three dimensions, for various values of mo: At 
relatively small times (i.e., larger than some microscopic scale but smaller than r m ) P c (t) decreases algebraically with 
the power 9q characteristic of the decay of the curve corresponding to mo = 0. As time t increases and exceeds r m , 
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however, one observes a crossover towards an algebraic decay with the power 9^ characteristic of the decay of the 
persistence probability for m = 1. (For comparison, in Fig. HJa) we also indicate the straight lines corresponding 
to algebraic decays with the powers 9 0lO c — to predicted in Sec. IIIIBI and reported in Fig. [2] and, c.f., Tab. Mil ) 
As the effect of a finite too is to introduce a time scale r m into the problem, it is natural to wonder whether these 
curves corresponding to different values of mo are characterized by dynamical scaling, i.e., if they collapse onto a 
single master curve after a proper rescaling of the time t and of the probability P c (t) involving r TO . The natural 
heuristic candidate for such a scaling form is P c (t) ~ r^ l 0o V(t/T m ), where, for consistency with the known behaviors 
for t <C r m and t 3> r m , one has V(x « 1) ~ x~ e ° and V(x ^> 1) ~ x~ 9 ^ . This scaling ansatz can be tested by 
plotting r^P c (t) or, equivalently, iru 6 " P c (t) , as a function of t/r m or, equivalently, m^/ K t where the exponents 
$o = 0o(d, d') and n(d) are the appropriate ones to the dimensionality of the model and of the manifold. In Ref. [ll| it 
was shown that the persistence probability of the global magnetization (d! = d) of the two-dimensional Ising-Glauber 
model obeys indeed such a scaling form with the expected value of k = 9[ s + f3/(vz) ~ 0.249. This numerical value 
follows form Eq. ©, where one uses the values of the exponents reported in Tab. [TT| together with the fact that 
in two spatial dimensions (3/v = rj/2. Similarly, we have checked (data not shown) that the same scaling behavior 
holds for the global magnetization of the three-dimensional Ising model, with the exp ected value of k — 0.34(1), 
which follows from Eq. © with 9[ s = 0.104(3) [Hj], P = 0.3267(10), v = 0.6301(8) [30], z = 2.04(2) 29]. However, 
the present numerical data indicate that this heuristic scaling ansatz does not capture the actual behavior of the 
persistence probability for the fluctuating magnetization of a manifold, neither in the case of a line in two dimensions 
nor of a plane in three dimensions. The data for a plane in three dimensions are shown in Fig. Eta), where we plot 
m -e (d=3,d'=2)/«p c ^ ag a function of m i/K t with = 3jd , = 2 ^ = 0(MC)( d = 3jd ' = 2) = 1.1 [see Fig. E|»] 

and k — 0.34 as derived above: clearly the curves corresponding to different values of the magnetization mo do 
not collapse onto a single master curve. We have carefully checked that the absence of data collapse is not due to 
finite-size effects. This lack of scaling can be qualitatively understood on the basis of the fact that the persistence 





10° F 




io- 1 ; 




io- 2 ; 




10 3 ■ 


Q_ 




10" 4 ■ 




10 -5 - 




10 6 




10 u 



10' 



10' 



1 3 1 4 1 o 5 





10 2 






ii 




T3 


10° 


CO 

II 








o 
CD 

O 


10" 2 


E 




O 

Q_ 


10" 4 



10 - ' 







ni =0.1 ' + 






m =0.3 






m =0.5 






m =1.0 □ 









10 - ' 



10 -2 10 _1 10° 10 1 10 2 10 3 



t m, 



1/k 



(a) 



(b) 



FIG. 4. Critical persistence probability P c (t) of the fluctuating magnetization of a plane within the three-dimensional Ising 
model for different values of the initial magnetization mo. In panel (a) P c is plotted as a function of time t and its long-time 
algebraic decay clearly shows a crossover, which is controlled by the initial value mo of the magnetization. The uppermost 
and lowermost data sets are the same as those presented in Fig. [3jb), whereas the associated straight lines indicate algebraic 
decays with the exponents Oq 1 ^ and 9^ s \ respectively, analytically estimated in Sec. IIII 51 (see Fig. [2] and, c.f., Tab. IIII| |. 
In panel (b) the data presented in (a) are scaled according to the heuristic scaling ansatz discussed in the main text: The 



9 (d=3,d'=2)/K 

g(MC) 



should collapse onto a single master curve when plotted as a function of t x ml/ K with 



data for P c (t) x m 

9o(d = 3, d' = 2) = ^q" 1 ^ (d = 3, d! = 2) = 1.1 and k = 0.34 [see Eq. P)]. However, this is clearly not the case and a significant 
dependence on mo is still observed in the resulting scaled curves. 



probability we have discussed so far is, in fact, a special case of a two-time quantity. Indeed, consider the persistence 
probability P c (t,t'), defined as the probability that the process X(r,t) does not change sign between the times t' 
and t > t'. In terms of this quantity, the persistence probability P c (t) studied above is given by P c {t) — P c (t, t^), 
where t^ is some non-universal microscopic time scale set, e.g., by the elementary moves of the dynamics. The 
scaling form of the correlation function (X(r,t)X(r,t r )) = C(t/T m ,t' /r m ) [see Eq. (TT71) ] implies straightforwardly that 

P c (t,t') = V(t/T m ,t' /r m ), where r m oc m^ 1 ^. In addition, the scaling behavior of the correlation function becomes 
independent of r m within the following two different regimes — which can be investigated analytically by taking, 
respectively, the limits r m — > oo and r m — > 0: (I) t' < t <C T m within which (X(v,t)X(r,t')) = C l (t/t') and (II) 
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T m <f'<i within which (X(r,t)X(r,t')) = C u (t/t'). Correspondingly, the persistence probability P c (t,t') takes two 
different scaling forms P: 

(I) P c (t,t / )=P I (t/t , )*~*' (t/t')- eoid ' d,) for t'<t<^T m , (39) 

(II) P c (t,t')=V n (t/t') t: ^' (t/t')- ^ d ' d,) for r m <t'<i. (40) 

In passing we mention that, interestingly enough, two analogous regimes emerge in the study of the persistence of 
fluctua ting interfaces, in which the role of r m is played by the equilibration time r eq ~ L 2 for a system of finite 
size L [3l|. As anticipated above, the case tuq = corresponds to regime (I) [Eq. (|3"9"|) ]. whereas m-o = 1 to regime 
(II) [Eq. (|40p]. While these two instances are well understood, our simulations for intermediate values of mo are 
likely to correspond to a third regime t' « r m « i which we are presently unable to describe analytically for generic 
D > because the correlation function of the process cannot be made stationary due to its residual dependence 
on T m . [Note that the introduction of the logarithmic time renders correlation functions stationary as long as they 
depend only on t' /t or, more generally, on ratios f(t')/f(t), see further below.] In particular there is no obvious 

reason why the scaling P c {t) ~ t~ 9a P(t/T m ) or, equivalently, P c (t) ~ mQ°^ d ^ K P c (m^ K t) — which well describes 
the behavior of P c (t) for D = [16| — should actually hold in general and therefore one should not be surprised by 
the fact that this scaling form turns out to be violated in the case of Fig- EJb) . Within the Gaussian approximation, 
however, one can heuristically understand why the case D = investigated in Ref. [l6| is indeed special in this 
respect. In fact, from Eq. (fTT|) one can see that only for D = the correlation function (X(r, t)X(r, t')) has the form 
(X(r,t)X(r,t')) = L(t'/T m )/L(t/r m ), with L(t) = t 1 / 2 ^ I (...,t), such that L(t) ~ t 1 / 2 for t < 1, whereas L(t) ~ t 2 
for i > 1 — we do not specify here the first argument of Jo as Id(x,u) is actually independent of x for D = 0, 
see Eq. (|18| . The persistence probability for such a process can be calculated exactly after a mapping onto a 
stationary process: P c (t,t') = (2/7r)arcsin[i(t'/T„ l )/L(t/r„ l )]. Therefore in the third regime t' -C r m < t one has 
P c (t,t') ~ (t' /Tm) 1 / 2 1 'L(t/r m ) which actually renders P c (t) = P c {t,t' = t„) ~ T m 9o P(t/ T m), i-e., the heuristic scaling 
form anticipated above, with the correct exponents 8 Q = 1/2 and 9^ =2 [la |. Finally, even though such a scaling form 
does not work for U > our numerical simulations indicate that data corresponding to different values of mo actually 
fall — to a rather good extent, see Fig. [5] — onto a single master curve obtained by plotting m °( d - 3 < d ~ 2 ^ k p c (t) 
as a function of m^ k t with 6*o(<i = 3, dl = 2) = 1.1 and k = 0.54 ^ n. However, the origin of this effective scaling is 
presently unclear and it surely deserves further investigations. 




FIG. 5. Effective scaling of the Monte Carlo data presented in Fig. |4ja) for the critical persistence probability P c {t) of the 
fluctuating magnetization of a plane within the three-dimensional Ising model. Data corresponding to different values of the 
initial magnetization mo collapse onto a single master curve when P c (t) x m -eo(d=3,d'=2)/fc ig plotted 

as a function of tm n 

with 6 {d = 3, d' = 2) = 1.1 and k = 0.54 / k. 



C. Comparison between analytical and numerical results 



The comparison between the analytical predictions of Sec. IIII Bl and the available numerical estimates of 6>o,oo 
discussed in Sees . IIV Al and HV Bl or reported in the literature is summarized in Fig. [5] and, for D — 1, also in Tab. IIIII 
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L> = 1 




d= 2 


d = 3 




0.78(5) 
0.64 


1.10(5) 
0.89 


ii(MC) 
foo 

/)("E) 
foe 


4.2(1) 
3.0 


3.7(1) 
3.5 



TABLE III. Monte Carlo and analytical estimates for the persistence exponents 8q and 9^ of the magnetization of a manifold 
of codimension D — 1 — i.e., a line in spatial dimension d = 2 or a plane in d = 3 — for an Ising model with Glauber dynamics, 
relaxing from an initial state with mo = and mo 7^ 0, respectively. The corresponding Monte Carlo estimates have been 

obtained as described in Sec. IIVI whereas the analytical estimates fl^oo = ^"Ooo x ^0,00 follow from the direct extrapolation to 
D — 1 of the results of Sec. IIII Bl summarized in Tab. [II] The values dj^j for d = 2 and d — 3 indicated here are reported, 
respectively, as gray and black vertical bars for D = 1 in Figs. [T]and[2] 



In particular, panels (a) and (b) of Fig. [5] compare the various estimates of 9^ and 9q, respectively, as functions 
of the codimension D, in spatial dimension d = 2 (grey lines and markers) and d = 3 (black lines and markers). 
As already pointed out in Refs. 0, [l^j , for vanishing codimension D = the agreement between the Monte Carlo 

estimate 0^q ^ (vertical bars in Fig. [2]) and the analytical estimates 6>^ S q (dash-dotted lines) of #oo,o is very good both 
in d = 2 and d = 3. For D = 1, instead, the comparison is less clear. On the one hand, our Monte Carlo estimate of 
9^ C \d = 3, d' = 2) = 3.7(1) in d = 3 (indicated in Fig. [2ja) as a vertical black bar — see Sec. HVBl for details and 

Tab. IIIip . is in a rather good agreement with the corresponding analytical estimate 9^ s \d = 3, d' = 2) = 3.5 reported 
in Tab. Mil [dash-dotted black line in Fig. Ufa)] and also with the result of the Gaussian approximation (solid and 

dashed lines). On the other hand, 9^ C \d = 2, d' = 1) = 4.2(1) in d = 2 (vertical gray bar in Fig. [U[a); see Sec. IIV Al 
for details) is significantly larger than the corresponding analytical estimate 9x ,(d = 2, d' = 1) = 3.0 reported in 
Tab. HI] [dash-dotted gray line in Fig. Ufa)]. In addition, the fact that 9 { ™ C) (d = 2,d' = 1) > 9 { ™ C) (d = 3, d! = 2) is 
at odd with what is observed for all the other numerical and analytical estimates presented in this work (c.f., Fig. [6] 
for the 0(n — > 00) universality class), i.e., that 9 OQ fi(d = 2,d' = 1) < 9 OQ: o(d = 3,d' = 2V As far as O for D = 1 
is concerned [see Fig. [3Jb)], we note that the preliminary Monte Carlo estimates of Ref. [l3| (black and grey circles, 
corresponding to 9^ AC \d = 3, d' = 2) ~ 0.88 and 9^ AC \d = 2,d' = 1) ~ 0.72, respectively) are essentially compatible 
with the values of 9^ gS> [dash-dotted line in Fig. [2jb); see also Tab. IIII| both in d = 2 (grey) and d — 3 (black). 

The Monte Carlo estimates #q MC - ) discussed in Sees. IIV Bl and IIV Al and reported in Tab. IIIII instead, turn out to 
be significantly larger than and marginally compatible with the previous ones in d = 3 and d = 2, respectively, and 
therefore they no longer agree with the corresponding analytical predictions provided by 9^^ . The reason for such a 
discrepancy might be that for D — 1 the actual contribution of the term of O(D) in fio x TZ^ 1 ^ [see Eq. (|3"8"]) ] is a rather 
large fraction of the O(D ) term (ca. 100% and 150% in d = 3 and d = 2, respectively), i.e., the dependence on D 
is so pronounced that a quantitatively reliable estimate of these exponents might require accounting for higher-order 
terms in the codimension expansion. 

We conclude this section by discussing briefly the case < ( < 1 (see Tab. |]|, that includes the instance of a line 
d' = 1 within the three-dimensional Ising model (i.e., D = 2) for which £ = ij/z ~ 0.016 (see Tab. ITTT) . For < £ < 1, 
our analysis indicates that the long-time behavior of P c {t) is independent of the actual value mo (e.g., mo = or 
mo = 1) and is characterized by a stretched exponential behavior (Tab.|I|. Unfortunately, in our simulations, we have 
not been able to observe incontrovertibly such a stretched-exponential law. This difficulty mirrors the one recently 
encountered in the numerical analysis of the persistence probability P c (t) of stationary processes characterized by 
two-time correlations with a power-law decay Cgt^i,^) ~ |£i — ^l - ^ [lH |33[. Although that context is rather 
different from the present, there it was observed that the convergence to the stretched exponential behavior — which 
is expected on the basis of the theorem by Newell and Rosenblatt 0] — is actually extremely slow. In addition, the 
pre-asymptotic behavior was shown to be increasingly important at the quantitative level as £ decreases [33[. Given 
these results and the extremely small value of C = 0.016 in the case of present interest it is not surprising that the 
stretched exponential is rather difficult to be observed. On the other hand, our numerical data indicate that the 
pre-asymptotic behavior of P c (t) is affected by the actual value of the initial magnetization mo- However, beyond this 
qualitative feature, we have not attempted a more detailed and quantitative characterization of this pre-asymptotic 
regime, which would require more extensive and dedicated simulations. 
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V. 0(n) MODEL IN THE LIMIT n -)■ oo 



Here we extend the previous analysis of the persistence probability of a manifold to the case of a n-component 
vector order parameter (p = (ipi, tp 2 , • • • , f n ) with n > 1 and model A dynamics [see Eqs. Q and ([5])]. Even though, 
strictly speaking, this model is not relevant for the description of actual physical systems, the fact that it can be 
solved exactly in the limit n — > oo provides insight into the effects of non-linear terms in the Langevin equation 
beyond perturbation theory. If the initial magnetization m does not vanish the original 0(n) symmetry of the model 
is explicitly broken and one has to distinguish between fluctuations Vv(x, t) which are parallel to the average order 
parameter m(t) — ((p(x,t)) and those i/j„{x.,t) which are transverse to it. The different fluctuations ip <7 (x,t) and 
V>7r( x iil are expected to have distinct persistence probabilities P%{t) and P£{t) which, on the basis of the results of 
Ref. [22| , can be shown to exhibit the same temporal crossovers as in Tab. U due to mo ^ 0. 



A. Persistence of the transverse modes in the 0(ri — > oo) model with finite mo 

We focus here on the limit n — > oo and we consider the transverse modes only, for which the response and correlation 
functions can be calculated exactly [22|. The wave- vector dependent response function is given by (see Eq. (68) in 
Ref. (H) 



^ (t>M) = ^Q,t') 



from which one obtains the correlation function (see Eq. (56) in Ref. [1^]) as 

C3(M') = (^(Q,i)^(-Q,* / )> = 2 / < dhR^t^R^t',^) 

Jo 

2(<t') £/4 



[(l + t/T m )(l + t'/T m )] 



_ e -Q 2 (* + I dt 1 (t 1 )- £ / 2 (l + t 1 /r m )e 2 « 2 * 1 



(42) 



where e = 4 — d. Due to the residual 0(n — 1) symmetry in the internal space of the transverse components, only 
correlation functions between the same component of the vector tp* can be non- vanishing. Accordingly, here and in 
the following we always refer to such correlations even if this is not explicitly indicated. For mo — > oo, r m ~ tn^ 1 — > 
and therefore one has 

i^(t>*',0=e- 0(t - t) (^) and 

V J t (43) 







Accordingly, the connected correlation function of the magnetization M 7r (r, t) of a manifold defined as in Eq. ([7]) with 
the substitution ip — > i/i* is 

<=i *" (44) 
=2c 2- D / 2 (tt') e/4 - 1 / 2 i /2 " (e+D)/2 C &uu x -^{y-u)- D l* , 



where f > t', y = (1 + x 1 )/2, and x = t' /t < 1 [cjj was defined after Eq. (|TB)) ]. The normalized process X^^i) = 
M 7r (r, t)/([M 7r (r, i)] 2 ) 1 / 2 is therefore characterized by the correlation function 

{X"(t, t)X*(r, t')) = (t/t')-( d+D V A A*(t'/t; d) , (45) 

which is of the form (fill) with -Foo (*/*') = ^ (*'/*! d), where 

A*(x; d) = x~ D ' 2j -^ - - , 46 
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(well defined for d > and < D < 2) and Hoo(d, d') = (d + D)/4. This expression for /ioo agrees with Eq. (|12l) in 
which 9 takes the value 9 = —/3/{vz) = — (d— 2 + r/)/(2z) characterizing the scaling form Eq. (J5J for the transverse 
modes (see, e.g., Refs. [H, [34j]) and which yields, in general, /i^ = 1 + [(d + -D)/2 — 2 + ?/]/z. For this specific model, 
?7 = and z = 2. According to Eqs. (f4"5"j) and (|46p the process X 1 is Markovian for vanishing codimension D — 0, with 
Moo = d/4. The corrections due to D ^ can be accounted for perturbatively by expanding A 7T (x; d) = x D ' 4 A 7r (x; d) 
for small D as 



A n (x;d) = 1 + DAl{x;d) +0{D 2 ), with 

Al{x-d) = --\nx-- C du «- w ' 2 h^ 
4 4 7 l-u 

After some algebra one obtains from Eq. ([2~5j) , 



(47) 
(48) 



y=(l+ a: -l)/2 



where 



8tt 



^L = l + DI d + 0(D 2 ) , 



-l+d/4 



dx- — 

o (1 



x d/2\3/2 



du u- 1+d ' 2 \n 



y 



l-u 



y={l+x- 1 )/2 



is a monotonically decreasing function of d, with Ii = 1.149 . . ., 1% = V2 — 1/2 = 0.9142 . . ., 7 3 = 0.867 . . . 
I4 = 1/n + 7r/4 — 1/4 = 0.853 . . .. For later convenience we indicate here also the values of Id for some d < 2. 
provides the estimate 



(49) 



(50) 

and 
This 



9W 

OO 



-[1 + D x I d + 0{D 2 )] 



(51) 



which, for D — 0, renders the value of the global persistence exponent reported for this model in Ref. [16j . As 
discussed in Sec. IIII1 a different estimate of 9^ can be obtained by expanding only the ratio IZd to first order in the 
codimension D, while keeping the full ^-dependence of the Markovian exponent fi^ = (d + D)/A, which corresponds 



to A n _= 1 in Eq. gSJ) and to A*{x; d) ^ A\(x; d) = Af(x; d) - (In x) /4 in Eq. gSJ. This yields 
Tl = n- D/d + 0{D 2 ), i.e., 



3 [2] 



9W = ^[l + Dx(I d -l/d) + 0(D 2 ) 



/loo x 72.: with 



(52) 



As expected, 9£i has the same small- D expansion as 6*Lo up to 0(D 2 ). 

For completeness and comparison we briefly present here also the expansion for 9q within the same 0(n — > 00) 
model, which was first discussed in Ref. [l3|. As the initial value mo of the magnetization vanishes, the original 0(n) 
symmetry is restored and there is no longer distinction between the transverse (tt) and longitudinal (a) fluctuation 
modes. The correlation function (X(r,t)X(r,t')) of the normalized process associated to the magnetization M(r,t) 



of the manifold can be obtained from Eq. (|44|) . with Cq = ^ q ^ (t, t') given by the limit for r m — > 00 of Eq. ((42 



CA(M') = 2(tt') e/4 e 



dti(ti)- 



e/2 e 2Q 2 t!_ 



(53) 



Interestingly enough, this expression of Cq(t, t') with r m = 00 [and, analogously, the expression for i?q(t, f'), see 
Eq. (|4ip] can be formally obtained from the one corresponding to r m = [see Eq. (|44p] by substituting in the latter 
e with e + 2, i.e., <i with d — 2. Accordingly, we can take advantage here of the results reported above for the case 
T m = in order to discuss the persistence properties of the manifold for mo = 0. In particular, for the normalized 
process X(r,t) = M(r,t)/([M(r,t)] 2 ) 1 / 2 , one finds 



{X{v,t)X{v,t')) = (X«(r,t)X*(r,t'))\ d ^ d _ 2 

= (t/t')- (d+D - 2 ^ 4 A 7T (f /t;d - 2), 



(54) 



on the basis of Eq. (gSJ). This result is of the form dTTJ> with /j, (d, d') = (d + D- 2)/4 and F (t/t') = A*(t'/t; d - 2). 
As expected, for d = 4 the equations above render the corresponding ones for the Gaussian model with mo = 0, 
which were discussed at the end of Sec. 1111 Al In addition, this expression for /iq agrees with the hyperscaling relation 
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(a) 




FIG. 6. Persistence exponents 8oo,o(d,d — D) of the transverse components of the global order parameter of a manifold of 
codimension D as a function of D, for d — 4 (uppermost solid and dashed curves) and d = 3 (lowermost solid and dashed 
curves), within the 0(n — > oo) model universality class with relaxational dynamics [Eq. Q]. Panels (a) and (b) refer to the 
case of non-vanishing and vanishing initial value mo of the order parameter, respectively, whereas the two different estimates 
@oo o [Eqs. (|5ip and (|55[) ] and 6^ [Eqs. (|52|l and (I56|l ] are indicated by the the solid and the dashed curves, respectively. For 
D > 2, i.e., C, > (shaded areas) the asymptotic decay of the persistence probability of this model is no longer algebraic. In 
panel (a), the uppermost solid and dashed curves (d = 4) are the same as those reported in Figs. [lib) and[2jb). The persistence 
exponent of the longitudinal component of the global order parameter for mo — is the same as the one for the transverse 
component [panel (b)], whereas for mo ^ and d = 4 (Gaussian approximation) it is given by the solid and dashed lines in 
Figs.[]Ja) and[2{a). The case of the longitudinal component for mo = was also discussed in Ref. [l3T ]. 



[analogous to Eq. (fl~3|) ] briefly mentioned after Eq. (|3 1 1) because r\ — and 9 = (4 — d)/A for the present 0(n — > oo) 
model [l||. On the basis of the mapping highlighted above, one can take advantage of Eqs. (|5"Tj) and (f5"2")) in order to 
determine the expansion of 6o(d, d — D) in the codimension D: 

[ o ] = + D x I d _ 2 + 0(D 2 )], (55) 

(which reproduces the expansion provided in Ref. [l3| for d = 4, whereas for d = 3 this gives a coefficient 2i/4 = 
0.287 . . . which corrects the value 0.183615 . . . reported therein for the first-order correction in D) and 

0[ 2 } = d-2 + D {1 + Dx [7 ^ 2 _ l/{d _ 2)] + 0{d2)} (56) 

In Fig.E^a) and (b) we report the estimates 9od for 9^, [Eqs. (J5T|) and ([52])] and ^ ' for 9 [Eqs. ((55} and (|56|) . see 
also Ref. [II]]) respectively, as functions of the codimension D. In the two panels 9^ and ^ are indicated as solid 
and dashed curves, respectively, for d = 4 (uppermost solid and dashed curves) and d — 3 (lowermost solid and dashed 
curves). Comparing the curves for d = 3 to those corresponding to d = 4 (Gaussian approximation), it turns out that 
the effect of non-Gaussian terms in the effective Hamitonian §5§ for n — > oo is to reduce the value of both 9oo and #o 
for generic values of D. This trend agrees with the one observed at the end of Sec. MI 51 both for the longitudinal and 
the transverse components of the order parameter and based on the results of a dimensional expansion around d = 4 
for the case D — and generic n. Even though a definitive statement in this respect would require a careful analysis, 
one might heuristically expect on the basis of the evidences collected here that, for a fixed codimension D ^ 0, the 
value of the persistence exponents #oo,o decreases as d decreases below the upper critical dimensionality d = 4 of the 
0(n —t oo) model, as it happens for D = 0. 



B. Relation between the spherical model and the 0(n — > oo) model 

We conclude this section by discussing the relation between the 0(n — > oo) model studied in Sec. IV Al and the 
spherical model (see, e.g., Ref. [35|). and its implications for the persistence probability. It is well-known that these two 
models are equivalent in equilibrium as there is a mapping between the corresponding free energies. This equivalence 
also extends to their equilibrium dynamics. However, when considering non-equilibrium properties some aspects of the 
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spherical model require a careful consideration. In particular, non-Gaussian fluctuations of the Lagrange multiplier 
which is introduced in order to impose the pure relaxational dynamics have to be accounted for when calculating some 
correlation functions of global quantities — which correspond to vanishing wave- vector Q = (see Refs. [13, 0|)- I n 
doing so, it turns out that the behavior of these quantities cannot be obtained as the limit for Q — > of the correlation 
functions for Q 7^ (local), given that the associated non-connected part (if non- vanishing) alters the scaling behavior 
compared to the case Q ^ 0. (Here and in what follows we assume that the system is spatially homogeneous.) 

Consider, for example, the correlation function of the local order parameter (i.e., of the magnetization). As long as 
the average of the local magnetization vanishes — which is the case if mo = — there are no differences between the 
connected and the non-connected correlation functions. However, as soon as the average of the local magnetization is 
non- zero (e.g., when the system is prepared in a magnetized initial state m ^ 0), the connected correlation function 
of the magnetization differs from the non-connected one only at Q = and the subtraction of the non-connected 
part alters its original scaling behavior. In turn, this subtraction affects also the scaling behavior of the response 
function. For a different observable, e.g., the energy, this difference in the scaling behavior might emerge also in the 
case of mo =0. Focusing here on the correlation and response functions of the order parameter of the spherical model 
(S), one finds that for mo = 0, the correlation and response function of global quantities can be simply obtained as 
the limit Q — > of the local quantities corresponding to a non-vanishing Q and that they coincide with the same 
quantities in the 0(n oo) model. For mo 7^ and local space integrals of the order parameter (i.e., Q 7^ 0, for 
which the non-connected part vanishes), instead, the correlation and response functions are given by (see Eqs. (2), 
(13) and (22) in Ref. 0) 

«5('>''.0 = (^)' /2 (^)' /2 e-"<-'», (57) 

where uj = ujq — 2 £™ =1 (1 - cosQ,), p = (4 - d)/2 = e/2 for 2 < d < 4, whereas p = for d > 4 (see Eq. (14) in 
Ref. EH] ), T m = cm,Q 2 (see Eq. (15) in Ref. [3]). The correlation function turns out to be, for t > t' 



2(U')P/ 2 

[(l+t/T m )(l+t'/T m )Y>' Z JO 



C^(t,t r ) =T c x ^-L r^e-^') / dt 1 (t 1 )-"(l+t 1 /r m )e 2 - tl . (58) 



In Eqs. (|57l) and (1581) . the superscript S indicates that the expression refers to the spherical model. Comparing 
Eqs. (|57|) and (|58|) with Eqs. (|4Tj) and (1421) . respectively, one concludes that, in the limit of small momenta (such that 
w = wq- Q 2 ), 

R^ (t,t') = R^(t,t') and C^ (t,t') = C^(t,t') , (59) 

up to non-universal factors, for all dimensions d > 2 (one can easily check that the equality holds also in the case d > 4, 
for which the Gaussian approximation is exact), and for all values of the initial magnetization m 7^ 0. We emphasize 
that the left-hand sides of Eq. (1591 refer to the order parameter of the spherical model, whereas the corresponding 
right-hand sides to the transverse fluctuations (tt) of the order parameter of the 0(n — > 00) model. These equalities 
justify the fact — already pointed out right after Eqs. (73) and (126) in Ref. [22[ and at the end of Sec. 5 in Ref. |24{ - 
that the asymptotic value of the fluctuation-dissipation ratio for local integrals of the order parameter in the spherical 
model is the same as the corresponding one for transverse fluctuations in the 0(n — > 00) model. 

Consider now the global order parameter, i.e., the integral over the whole space of the local order parameter, 
corresponding to the Fourier mode with Q = 0: its statistical average does not vanish for mo 7^ and therefore its 
connected two-time correlation function differs from the non-connected one and cannot be simply obtained as the 
limit for Q — > of the same quantity for Q ^ 0. The explicit expression for the correlation Cq = o and response 
Rq=o functions of the global order parameter for d > 4 and arbitrary value of the magnetization mo were reported 
in Ref. [see Eqs. (48) and (47) therein]: 



R S Q=0 (t>t',t')=( 1 J ^l^) 3/2 , (60) 



and (written in a slightly different form compared to Ref. (24|) 

" l(l + t/T,„)(l + t'/T m )f' 2 

By comparing these expressions with Eqs. ([57)1 and (|58|). one realizes that — as anticipated above — the former are 
not given by the limit Q — > (i.e., ui — > 0) of the latter. On the other hand, the expressions for Rq =0 and Cq =0 
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for the spherical model in d > 4 are the same (up to an irrelevant multiplicative factor for C) as the one for the 
longitudinal fluctuations (a) with Q = of the order parameter of the 0(n — > oo) model, given by Eqs. (58) and (59) 
of Ref. [H: 



RS(t, t') = R°(t, t') and C* (t, t') = t') 



(d>4) 



(62) 



where, as the superscripts indicate, the left-hand sides refer to the spherical model whereas the right-hand sides to 
the longitudinal fluctuations in the 0(n —> oo) model. (Note that for the case d > 4 presently discussed there is no 
difference between the longitudinal fluctuations of the 0{n) and 0(1) model, and therefore the expressions for i?q =0 
and Cq =0 can be read from Ref. in which the latter, i.e., the Ising model is actually studied.) Even though we 
have proven the relation (|6"2"|) only in the case d > 4, one heuristically expects it to be valid also for d < 4, as it 
was the case for the analogous relation (|5"9"|) . The expressions for the correlation and response function of the global 
magnetization of the spherical model for 2 < d < 4 and a generic value of the initial magnetization mo (i.e., r m ) can 
be found in Ref. (24|. While the expression for -Rq =0 is explicitly given by Eq. (58) therein, 



Rg =0 (t>t',t') 



l-d/4 



1 + t'/Tr, 
1 + t/T m 



1/2 



t/r„ 



1 + t/r n 



t ,\d/2-l 



(63) 



the corresponding expression for the correlation function Cq =0 (i,i') is significantly more involved and it has been 
worked out explicitly only in the asymptotic regime t' <C t. For the 0{n) model, instead, the response and correlation 
functions RQ= a {t,t') and Cq =0 (i,i'), respectively, have been calculated for d < 4 only at the first order in the e- 
expansion around d = 4 (with e = 4 — d) and in the limit of large magnetization, i.e., for r m —¥ 0. These expressions, 
reported in Ref. [22[ can be compared in the limit n — > oo with the first term of the expansion around d — 4 — e of 
the corresponding result for the spherical model. In particular, focusing on the response function for n — > oo, one has 
(see Eqs. (88) and (89) in Ref. |l|) 




0(e 2 



(64) 



which is indeed equal to the first-order expansion of Eq. (|63l) around d — 4 and for r m — 0. In order to do an analogous 
comparison for the correlation functions Cq_ (i, i') and Cq_ (t,t') for d < 4 one would have to calculate the limit 
n — > oo of the results for Cq_ (t, t') presented in Ref. [22] for r m = [in particular, see Eqs. (93), (B31), (B24), and 
(B29) therein] and compare it with the limit r m — > of the correlation function Cq_ (i, t') of the spherical model 
explicitly presented in Ref. jUJ for the case t' <C t. However, this somewhat lengthy calculation can be avoided by 
noting that this equality up to C(e 2 ) is implied by the one between the response functions of the spherical and of the 
0(n — > oo) model up to C(e 2 ), together with the fact that the asymptotic value of the (fluctuation-dissipation) ratio 
[d t 'C(t,t')\/ R(t,t') for t > t' and r rn -> is the same in the two models [see right after Eq. (105) in Ref. [jj|]. 

Summing up, the results of Refs. 0, [22| - [24| suggest the equalities Eqs. (I5T)|) and (|6"2")l for the order parameter 
response and correlation functions of the spherical and of the 0(n oo) models relaxing from an initial state with 
mo 0. Whereas Eq. (|59[) can be explicitly checked for d > 2 and all values of the initial magnetization mo 7^ 0, 
Eq. (|6"2")l is proven for generic values of r m only for d > 4, while in the case T m = also for 2 < d < 4 but only in the 
limit t' <C t. If, as it is likely, this relation extends to the remaining cases, it would be interesting to understand its 
deeper motivation as it connects, together with Eq. ([5^| . different degrees of freedom in different models. 

According to the correspondence highlighted above, for a spatially constant non-vanishing initial value of the order 
parameter mo, the persistence properties of the global order parameter of the spherical model are (up to non- universal 
factors) the same as the one of the global longitudinal fluctuations (cr) of the order parameter in the 0{n — > 00) model. 
However, the persistence properties of the global order parameter of a manifold of non- vanishing codimension D =/= 
in the spherical model are (up to non-universal factors) the same as the ones of the global transverse fluctuations (tt) 
of the order parameter of the same manifold in the 0(n — > 00) model. A consequence of this correspondence is that 
the small codimension expansion which typically allows for a perturbative access to the persistence exponent of the 
manifold (D ^ 0) within the spherical model is not an expansion about the persistence exponent of the global order 
parameter of the entire system (D = 0), as highlighted by the fact that the former and the latter actually involve 
different degrees of freedom (ir and cr, respectively) in the corresponding 0(n 00) model. 

In passing we mention that the result reported in Ref. for the global persistence exponent #oo of the spherical 
model (i.e., with D = and mo 7^ 0) is incorrectly based on the expression of the (connected) correlation function 
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C m (t, t') of the global magnetization for t ^> t' reported in Eq. (8.108) of Ref. (23|. Indeed, the knowledge of C m (t ^> 
t' , t') gives access only to the associated exponent — d/4+1 within the Markovian approximation (which is correctly 
reported in Ref. [Hj]), whereas 9^ 7^ fi^ is actually determined by the full functional form of C m (t, t') for generic values 
of t and t' , which is not explicitly provided in Ref. (23[. The correspondence discussed above conveniently provides 
some information on the global (i.e., D = 0) 9^ for spherical model with mo ^ on the basis of the corresponding 
results for the 0(n — > 00) model presented in Ref. [l6j]: Indeed, it implies that 9^ — 9^, in which the left-hand side 
refers to the spherical model whereas the right-hand side to the longitudinal fluctuations of the 0(n — > 00) model. For 
the 0(n) model one finds 9 a x = ^xW, where ^ = l+d/(2z) &ndTZ a = l + e[(0.115 . . . + n0.131 . . .)/(8 + n)]+C(e 2 ) 
with e = 4 - d 0, so that, in the limit n4 00, /x^ -> 1 + d/4 (z 2 (H) and TZ a ->■ 1 + e x 0.131 . . . + C(e 2 ). For 
the spherical model this implies that 9^ = (1 + d/4) x [1 + e x 0.131 . . . + 0(e 2 )] which indeed shows a non-Markovian 
correction. Even though it would be nice to have a direct check of this prediction on the basis of the results of Ref. [23[ 
for generic d, the direct calculation of 9^ turns out to be rather involved and it is beyond the scopes of the present 
study. 



VI. CONCLUSIONS AND PERSPECTIVES 



In summary, we have investigated both analytically and numerically the persistence probability P c {t) of a d'- 
dimensional manifold within a d-dimensional system which relaxes at the critical point from an initial state with 
non- vanishing value of the magnetization or, more generally, order parameter mo- Such a persistence probability 
is defined as the probability that the fluctuating order parameter does not cross its average value up to time t. 
Depending on the value of the parameter Q = (D — 2 + rf)/z, with D — d — dl we found that, as in the case m = 0, 
the long-time decay of P c (t) is (i) exponential for £ > 1, (ii) stretched exponential for < C < 1 and (hi) algebraic 
for £ < 0. While in the first two cases (i) and (ii) the asymptotic behavior of P c (t) is not affected by a finite value 
of mo, in the third case £ < we demonstrated that P c (t) exhibits a temporal crossover between an early-time 
and a distinct late-time algebraic decay, which are characterized by two different exponents 9o(d 1 d l ) and 9 00 (d,d'), 
respectively, with 9 00 (d,d') > 9o(d,d') (see Tab. H|. Analogously to the case D = 0, the crossover is controlled by 
the time scale r m oc m l ^ K . The analytic determination of the associated exponents #0,00 (d, d') is rather non-trivial 
already within the Gaussian approximation (which becomes exact for d > 4) because the stochastic process under 
study turns out to be non-Markovian for D 7^ 0. In order to calculate 9oo(d, d') (9o was already studied in Ref. [l3j]) 
we performed a perturbative expansion up to the first order in the codimension D = d — d! of the manifold [see 
Eqs. (f2l)|). (|2"8")l , (1331) and (|3~4"1) ]. Then we presented a perturbative approach which allows one to calculate the effects of 
non-Gaussian fluctuations in a dimensional expansion around the space dimensionality d — 4. Combining these two 
expansions we obtained analytic estimates of 9oo(d,d') and 9o(d,d') up to order 0(D,e) [see Sec. IIII 51 and Fig. [2], 
In order to assess the reliability of these analytic estimates and to go beyond the perturbation theory, we studied the 
critical relaxation of the Ising model with Glauber dynamics on a e?-dimensional hypercubic lattice. In particular we 
computed the persistence probability of the magnetization of a line in d = 2 and of a plane in d = 3 — corresponding 
to codimension D = 1. In both these cases we observed a temporal crossover in P c (t) between two algebraic decays, 
as predicted by our analytical investigation and we determined the numerical estimates of the associated exponents 
9^ C \d = 2,d' = 1) and 9^ C \d = 3,d' = 2) as well as more accurate estimates of 9^ C \d = 2,d' — 1) and 

9 < ^ AG \d = 3, d' = 2). In the case of vanishing codimension D = — primarily investigated in Ref. [l6| — the 
agreement between the corresponding numerical and analytical estimates of #o.oc are very good both in two and 

three dimensions, as summarized here in Fig. [U(b). For D = 1, instead, 9^ C \d — 3, d' — 2) is in rather good 

agreement with the corresponding analytical estimate, whereas 9^ C \d = 2,d' — 1) is significantly larger than the 
value predicted by our perturbative calculation — see Fig.^a) and Sec lIVCI for the comparison. In addition, it turns 
out that 9^ c (d = 3, d! = 2) < 9^, c (d = 2, d' = 1) while all the analytical approaches presented here (including the 
analysis of the 0(n —> 00) model) suggest that the converse should be true. Our numerical simulations also unveiled 
a non-trivial scaling of the persistence probability with the characteristic time r m , which remains to be understood. 
These two latter intriguing features of the persistence probability of the manifold surely deserve further investigations 
beyond the preliminary one presented here. 

Finally, we have complemented our analysis of the persistence properties by a thorough comparison between the 
non-equilibrium dynamics of the 0(n — > 00) model and of the spherical model. While they are strictly equivalent as 
far as their equilibrium properties are concerned, we have pointed out that such an equivalence has to be carefully 
qualified when discussing their non-equilibrium dynamics. In particular, in the case of the critical relaxation from an 
initial state with non- vanishing order parameter, an unexpected connection emerges between the local order parameter 
of the spherical model and the transverse components of the order parameter in the 0(n — > 00) model as well as 
between the global order parameter of the former and the longitudinal components of the latter — see Sec. IV Bl for 
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details. This connection is related to the fact — pointed out in Refs. [23|, l24[ — that within the spherical model the 
correlation function of global quantities (corresponding to a vanishing wave- vector Q = 0) can not be obtained as the 
limit Q — ► of the associated local correlations (corresponding to Q ^ 0) . 

In view of the results presented here, it would be certainly interesting to analyze the consequences of a finite initial 
magnetization mo on other relevant properties which characterize the temporal evolution of the thermal fluctuations 
Sm(t) of the magnetization of a manifold. As an example, it was recently shown in Refs. [13, HH that the longest 
excursion l max _(t) between two successive zeros of a stochastic process up to time t is an interesting quantity which 
characterizes the "history" of the stochastic process and the asymptotic behavior of which depends qualitatively on 
the value of the persistence exponent of the process being larger or smaller than a certain critical value 9 C (3tJ • Given 
that, for the systems studied here, 0o(d,d') < 9 C < oo {d,d!) one expects an intriguing dynamical crossover in the 
growth of the average (Z max (i)) [37[, which certainly deserves further investigations. 
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